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Abstract. We give a non-iterative solution to a particular case of the four- 
point three- views pose problem when three camera centers are coUinear. Using 
the well-known Cayley representation of orthogonal matrices, we derive from 
the epipolar constraints a system of three polynomial equations in three vari- 
ables. The eliminant of that system is a multiple of a 36th degree univariate 
polynomial. The true (unique) solution to the problem can be expressed in 
terms of one of real roots of that polynomial. Experiments on synthetic data 
confirm that our method is robust enough even in case of planar configurations. 



1. Introduction 

We first recall some definitions from multiview geometry and formulate the prob- 
lem in question, see [3l[8] for details. A pinhole camera is a triple (O, n,p), where 
n is an image plane, p is a central projection of points in 3-dimensional Euclidean 
space onto 11, and O is a camera center (center of the projection p). The focal 
length is the distance between O and H, the orthogonal projection of O onto 11 is 
called the principal point. A pinhole camera is called calibrated if all its intrinsic 
parameters (such as focal length and principal point's coordinates) are known. 

Problem 1 (The four-point three- views pose problem). Let us consider three cal- 
ibrated pinhole cameras with centers Oi, O2, O3 and four scene points Pi, ... ^ Pi 
being in front of the cameras in 3-dimensional Euclidean space, see Figure [TJ In 
every local coordinate system Ojxyz, associated with the jth camera, the homo- 
geneous coordinates {xji yji l)"^ of all P^'s are only known. One must find the 
relative position and orientation of the second and third cameras with respect to 
the first one and the coordinates of all P^'s in Oixyz. 

Remark 1. It is well-known that the solution to Problem [T] is only defined up to an 
overall scale. That is, multiplying the coordinates of all the reconstructed points by 
any A > 0, we get another solution to Problem [TJ In order to resolve this ambiguity, 
we assume from now on that the length of the baseline joining Oi and O2 is given 
and equals d. 

Problem[T]is known to be slightly over-determined and in general to have a unique 
solution [31[TU] . The only iterative approach to the problem has been proposed in [5] . 

In the paper presented, we consider a light (but still highly nontrivial) version 
of Problem [T] when the three camera centers Oi, O2 and O3 are coUinear. We use 
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Figure 1 . To formulation of the four-point three- views pose problem 

the method of paper [7] to derive a relatively simple system of three polynomials 
in three variables. It is shown that the eliminant of that system is a multiple of 
a 36th degree irreducible univariate polynomial. It is conjectured that in general 
one of its real roots encodes a solution to Problem [T] Experiments show that our 
approach is comparable in accuracy with the existing five-point solvers. 
Here we outline our algorithm: 

(1) Transform the initial data Xji, yji in such a way that Xji = yji = Xj2 = 0, 
j — 1, 2, 3, applying rotations ([T]). 

(2) Derive from the epipolar constraints six polynomials by and one 
polynomial /™'^ by p3|) . 

(3) Construct system of three polynomials in three variables. 

(4) Using resultants and polynomial divisions, derive a 36th degree univariate 
polynomial S. 

(5) For every real root of S, find the matrices i?*-^' , i?'-'^' and [t] x by (O and JTU]) 
respectively. 

(6) Find the coordinates of O2, O3 and Pi by ((34l) . p5|) and (|36p respectively. 

The paper is organized as follows. In Section [21 we describe in detail our algo- 
rithm. In Section [3l we make a comparison of our algorithm with the five-point 
Li-Hartley solver [6] on two sets of synthetic data. In Section |4j we discuss the 
results of the paper. 



2.1. Transformation of the initial data. Initial data for our algorithm are the 
values Xji, yji for j = 1,2,3, and i — 1,...,4, as well as the distance d (see 
Remark [T]) . 

Without loss of generality, we can set Xji — yji — Xj2 ~ for j = 1, 2, 3. If Xji, 
yji and Xj2 differ from zero, we rotate Ojxyz to Ojxyz with the matrix 



2. Description of the algorithm 



Pj — Rip.Re-R^, 



(1) 
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where R^p■ , Re^ , Rip^ are rotations through the angles Lpj , 6j , respectively around 
the axes OjZ, Ojy, OjX respectively. The angles can be expressed by the formulas: 



tan lpj = 
tan 6i - 



Xji 



tan ipj 



{yj2 - V3i)\l^ + x]i + yn 



One then verifies that Xji — yji — Xj2 — 0. Besides, the depth of Pi from the 



camera center Oj is 



1 + + yl > 0, 



i.e. the point Pi remains to be in front of all the cameras. 

We will see that the above transformation of the initial data, being quite simple, 
noticeably simplify our further computations. 

2.2. Epipolar constraints and essential matrices. From now on we assume 
that the three camera centers Oi, O2 and O3 are collinear. 

Denote by (cc^j y'^^ l) and t = (tx ty t^) the directing vectors of OjPi 
and O1O2 respectively in the world coordinate system Oixyz. It follows that the 
directing vector of O1O3 is at, where a is either +1 or —1. 

Let Cartesian coordinates of Pi and Oj be {x p. , yp. , zp. ) and (xq ,yo t ) 
respectively. Then, it is easy to see that (j = 2, 3; i = 1, . . . , 4) 



{Xli - x'ji) tzZp^ = {tx - tzx'ji) ZOj , 

{yii - y'ji) tzzp^ = {ty - t^y'ji) zo^ ■ 

Eliminating zp. , Z02 and from this system, we get 

tx{y'2i - yii) + ty{xu - x'^i) + tz{yi,x'2, - xny'^i) = 0, 
txiv'^i - yii) + ty{xu - x'^i) + tz{yux'^^ - xuy'^i) = 0. 

Consider matrices i?*^^^ G S0(3), j — 2,3, rotating OjXyz so that 



(3) 



(4) 





R''' y.. ■ (5) 



where ujji is called the depth of the point Pi from the camera center Oj . Then, it 
is easy to see that Q is equivalent to the so-called epipolar constraints: 



xu yu 1) E^-^^ \y2^\ = 0, (6) 




where i?^^-* ~ [^Jx^^-''' is called the essential matrix and [t]x is the skew-symmetric 
cross-product operator. 
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2.3. Seven polynomials in six variables. Our approach is based on the follow- 
ing well-known result. 

Theorem 1 ( 1 ). If a matrix R G SO (3) is not a rotation through the angle tt 
about certain axis, then R can be represented as 

^ /l - + w^(l - w^) 2{vw'^-~u) 2w{u + v) \ 

— 2{vw'^+u) -w'^{l-v'^) 2w{uv-l) , (7) 

\ 2w{u~v) 2w{uv + l) l + i? -w'^{l + v'^)) 

where A = 1 + + w'^{\ + v"^). 

Let Oxyz and OXY Z be two coordinate systems with a common origin and of 
the same handedness. Then the Euler angles, transforming Oxyz to OXY Z , are 
defined as [5]: 

• if is the angle between the x-axis and the line of nodes, i.e. the line of 
intersection of the xy and the XY coordinate planes. 

• 9 is the angle between the z-axis and the Z-axis. 

• -0 is the angle between the line of nodes and the X-axis. 

Remark 2. By direct calculation one verifies that the Euler angles ip, 6 and %p can 
be expressed in terms of the parameters u, v and w as follows: 

u + V 

ip = arctan ■ 



1 ~ uv 

U — V 

lb — arctan , (o,) 



±2 arctan.. + 



Lemma 1. Let ipj, 9j and ■0j be Euler angles parameterizing the matrix i?^-'-'. In 
case of collinear camera centers Oi , O2 and O3 , the following identity holds: 

tan p2 = tan 953. 

Proof. Due to the condition Xji = yji = 0, j = 1,2,3, the z-axes of all Ojxyz 
intersect in the only point Pi. Since Oi, O2 and O3 are collinear, the lines of nodes 
of 02xyz and O^xyz are identical and hence, by definition, ip2 = p^ mod tt. □ 

Lemma 2. Let E'^'^\u,v,w) be an essential matrix subject to system If 

«; = --7'' (9) 

Vj yi2Wj{vj + Uj) + yj2Wj{vj - Uj) - 2yi2yj2Uj 



u 



] yi2{vj + Uj) - yj2{vj - Uj) + 2yi2yj2VjWj 



then £:(J) (u^. , v'^ ,10'^ ) = -E^^^ (u^ , ,w^). 

Proof. Let E'-^^ ~ [ijxi?'-'' be subject to Then, from the first two equations 
of ® (for i = 1,2) we find 



, R^^^(R^^^y,2 + R^^^)yi2 



{^23^12 - ^1^^22)1132 + Rii{Rz2yj2 + R'ii)yi2 

R[y + Vj 



(10) 
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where R^j} is the (i,fc)th entry of R^^\ Substituting these values into the ma- 
trix [t]x, we get = (uj , Vj, Wj). By a straightforward computation, we find 
that the equation E^-^^ {uj , Vj , Wj ) -fi?'-'^ {u'j , v'^ , ) = has the only solution ^ . □ 

Since system Q is linear and homogeneous in t^, ty and tz, it has a nontrivial 



(11) 



solution if and only if rank [ q{3) ] 2, where 



qU) = " yi2 -x'j2 yi2x',2 

y'js - yi3 - ^'j3 yi3x'j3 - xi^y'^^ 

Let i = i[a,0,-t) number triples from {1,2,3,4} so that j(i,2,3) — Ij *(i,2,4) — 2, 
*(i,3,4) — 3- Then, it follows that detQp^ — 0, where Qp^ is a 3 x 3 submatrix 
of Q*--'-' corresponding to the rows a, /3 and 7. Thus, we get a system 

#=detQF) A^^,.^./.^i. ^P^ ,=2,3, z = 1,2,3. (12) 

where Aj = 1 + + li'j (1 + and, for instance, LOja is the depth of the point Pa 

from Oj (cf. (O). One verifies that in general each f^''^ is an irreducible polynomial 
in Wj, Vj, uij of the sixth total degree. 

Remark 3. In (|12p . we assume that Wj 7^ 0. No solution is lost here, since Wj — 
means that the matrix i?'^' is a rotation about the 2;-axis, which may only occur if 

In addition to (|12p. we can also obtain more polynomial equations by treating 



the rows of both Q*^^) and Q^"^-*. For example, consider a submatrix Q™"' of 
consisting of the first two rows of Q*^^'' and the second row 

ofQ(3)^ Then, 



lg(3) 



= det Aj A3 0.21(^22^32 ^ (^3^ 

^21 

where we assume that 0:21 ^ 0. In general, f™^^ is an irreducible polynomial in U2, 
V2, W2, M3, ^^3, W3 of the seventh total degree. 

2.4. Two polynomials in three variables. In this subsection, we deal with the 
polynomials fj^''^ defined in (fT2|). We are going to eliminate the variables W2 and 
W3 from them and then simplify the obtained polynomials using the transforma- 
tions ©. 

Notice that every polynomial /p'' is of the second degree in Wj, i.e. 

/(^■) =ap'u;2+5Wy,^.+c|^'\ (14) 
where a'f\ b\^\ c^"* ^^'^ polynomials in Uj and Vj. Consider a matrix 



a^^ hf cf 
\af hf cf 
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Then, system (fT2|) (for each j) has a solution if and only if det F^-'^ = 0, i.e. we get 
two polynomial equations: 

g<^i){uj,Vj) = detF(j') =0, J = 2,3. (15) 

The polynomials g*^^-' and g^^^ are of the 10th total degree in the variables U2, V2 
and U3, V3 respectively. We can reduce the total number of variables to three by 
introducing a new variable 

1 - U2V2 1 - ^3^3 

where the second equality holds due to Lemma [T] Note that the variable s is 
unchanged under the transformations (|9]), i.e. 

_ U'j + v'j _ U, +v, _ 



Substituting 
to P^ . we get 



1 - U'.V'^ 1 - UjVj 



V, = (16) 



By Lemmadl h^J^{u'j,s') ~ h^^'>{—Uj ^,s) = 0. Moreover, /i*^^^ has a special sym- 
metric form: 

12 



/^0)=.f ^pO-)[,12-.^(_^^.).-12]^ (17) 



k=0 



where are 6th degree polynomials in s. Due to the above symmetry, we can 



introduce a new variable Uj — uj — Uj , and then transform hj to the polynomial 

u=) = j2p^^^, (18) 

k=0 

where 

6 

Ai) _ \ " „(i) J,. , 

■f^2fc ~ / ,-P12-2i^».fc' 

EU) 
Pll-2i^i+l/2,fc+l/2, 



P2k+l 



i—k 

and Ki^k denotes the sum of two binomial coefficients: 

As a result, we have two irreducible polynomials h^^^ {u2, s) and /i^^^ (£13, s) of the 
12th total degree each. In order to solve the problem, we need at least one more 
polynomial in the variables U2, us and s. 



A SOLUTION TO THE 4P3V POSE PROBLEM IN CASE OF COLLINEAR CAMERAS 7 



2.5. One more polynomial in three variables. Let us consider the polynomial 
jmix (jggQg(j (jT3|) . As in the previous subsection, we are going to eliminate 
the variables W2 and W3 from J™"' and then simplify the result using transforma- 
tions dg. 

First, we notice that is of the second degree in the variable ■W2'- 

^mix ^ „mix^2 _^ ^mix^^ ^ ^mix^ ^21) 

where a™"', 6™"^, c™"' are polynomials in the remaining five variables. Consider a 
matrix 







6f) 




^mix ^ 




.(2) 
"2 






Umix 


^mix 





The three polynomials f['^\ f^'^ and y™'^ have a common solution if and only if 
det F™''' = 0. This yields 

det 

/™« ^ = (22) 

(2/12 - y22)w2 + 2(y?2 + 2^12^12 + 2/|2)"2t^2 + (y?2 " y22)^'2 

One verifies that the polynomial f"^^^ is in turn of the second degree in W3, i.e. 

jmix ^ ^mix^2 ^ gmix^^ ^mix^ (23) 

where a™™, 6™"^, c™"^ are polynomials in U2^ i'2i '^a and W3. Similarly, we define 
the matrix 

pmix _ (3) , (3) (3) 
-1 flj ^2 ''2 ' 

^^mix ^mix ^mix^ 

and from the constraint det F™^^ — obtain the polynomial 

det F™''^ 

mix ^'^^ ^ Q f24) 

~ (y?2 - yi2)"3 + 2(y?2 + 2y?2y|2 + y|2)"3t'3 + {yh ~ vl2>l ^ 

Substituting ^B)) to 5™"^ = 5™"'(m2, W2, W3, ws), we get 

(1 + su2)^(l + su^if g'^'''{u2, TT^, "3: -rr^) 

n (u2,M3,s)- — — — 2vrr"; — 2T7 ; tt ; ^ — ■ i^s) 

(1 + u^){l + u%)[xi'i + yi3s)(a;i4 + yus) 

Remark 4. The denominators in (P^ . and (^5]) are supposed to be nonzero. 
By a straightforward computation, one verifies that in general their roots do not 
give a solution to Problem [TJ 

The polynomial ft,™"' has a special symmetric form: 
4 4 

ft-'- = uiut ["2"' + (-"2)'=-^][^^r" + (-K3)'"~^], (26) 

fe=0 in=0 

where are polynomials in s. Substituting Uj = Uj — u~^, we transform ft™"' 
to the polynomial 

2 2 

^"'^ = EEpr™^2<, (27) 

fe=0 m=0 
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where 



2 



-mix 
P2k.2r, 



/ , P4-2i,4-2j^i,k>Cj,m: 
i=k j=m 
2 1 

-mix \ ^ \ ^ ^mix 

P2k,2rn+1 ^ 2^ 2^ P'l-2i,3-2j ^i.k + 1 /2,m+l /2, 
i=k j=m 
1 2 

P2k+l,2rn = 2^ 2^ P3-2i,4-2j^i+l/2,k+l/2>^j,m, 
i—k j—rii 
1 1 

P2A;+l,2m+l = / . j^3-2z,3-2.?^i+l/2,A;+l/2>^j + l/2.m+l/2i 

i—k j—rn 



(28) 



and >Ci^k is given by ((20)) . 

2.6. Some notions from elimination theory. In this subsection, we briefly 
recall some notions from elimination theory, see [2] for details. 

Given an ideal J — (/i, . . . , fm) C C[xi, . . . , the Ith elimination ideal Ji is 
defined by J; = JnC[x;+i, . . . , x„], i.e. J; consists of all consequences of fi — ■ • ■ — 
fm = eliminating the variables xi, . . . , xi. The generator of Jn-i is called the 
eliminant of J in the variable a;„ . 

As is well-known, the best way for finding elimination ideals is computing a 
Grobner basis of J with respect to the pure lexicographic ordering xi > ■ ■ ■ > 
Xn- However, in many cases this computation is practically impossible because 
of bounded computer resources. For these cases some roundabout ways, such as 
theory of resultants, should be applied. 

p q 

Let /i = X] s-nd /2 = ^i^i-i where Oi, bi £ C[x2, . . . , x„], be polynomials 

4=0 i=0 

in C[a;2, . . . The Sylvester matrix Syl^.^ (/i, /2) = (sfci)01i /i ^^'^ h 

with respect to xi is defined by 

Sp+q-j + l,j+i — bq^i, 

and all other entries of Syl^.^ (/i, /2) are equal to zero. The resultant of /i and 
/2 with respect to xi, denoted 7^m(./i,/2), is the determinant of the Sylvester 
matrix, i.e. TZx^ (/i , /2) — det Syl^.^^ (/i , /2). The following lemma states the relation 
between resultants and elimination ideals. 

Lemma 3 ([2]). Let /i,/2 G C[a;i, . . . , a;„] have positive degree in xi. Then, 

7^.,(/l,/2)G Ji. 

2.7. Thirty-sixth degree univariate polynomial. At this step, we have derived 
the following polynomial system in the variables U2, and s: 

hC^^ = [4]s^ul + [5]sul + [6]ui + [6]ul + [6]ul + [6]u2 + [6], 

^ [A]s^ul + [5]sul + [6]ut + [6]ul + [6]ul + [6]u3 + [6], 

/j""'^ ^ [3]sulul + [A]ulu3 + [4]u2ul + [4]u^ + [4]u2U3 + [Muj + [4]u2 + [4]% + [4], 

(29) 

where [n] means an nth degree polynomial in the variable s. 
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Let US consider an ideal J — {h^^\ h^^\ h™^^) C C[u2, "is, s] and denote its elim- 
ination ideals by Ji = J n C[u3, s] and J2 = J fl C[s]. 
By direct computation, we find 

TZu, {h^^^ , h'^n = r(s2 + if Sl e Ju (30) 

where r is a 28th total degree polynomial in the variables U3 and s, 52 and S3 (see 
below formula (|32)) ) are 4th degree polynomials in the variable s with the known 
(rather cumbersome) coefficients. For example, the trailing coefficient of Sk is 

j/fc2 [xi4Xk3{yi2 - yiz){yk2 - Vki) - xi3Xk4{yi2 - yi4){yk2 - yks)] 
X [a;i3(yi2 - 2/i4)(a;fe4 - Xk3Xk4 - ykzykA + yli) 

+a;i4(yi2 - yi3){xl3 - XksXkA - ykzykA + yfcs)] • (31) 

Further, 

nu,{U^\nu,{U^\h^'^))=S{s^ + iy^s\x^3+yi^s)\x,i + y^^sf^ e J2, 

(32) 

where 5 is a 36th degree irreducible polynomial in s. 

Lemma 4. The eliminant of the ideal J in the variable s is a multiple of the 
polynomial S, i.e. J2 = (Sp), where p £ C[s]. 

Proof. Let sq be a root of S. Since sq is not in general a root of the leading 
coefficients of /i^^^ and h'^^^ (denoted by [4]s^ in (HH)), it follows from the Extension 
Theorem [2] that there exist such U2fl and U3.0 that ("02,0: "3,0, so) is a solution 
to system ([29| . This exactly means that sq is a root of the eliminant of J in the 
variable s. □ 

By Lemma 21 a root of the polynomial S can be extended to a solution of 
system (|29)) . Moreover, we conjecture that in general every solution of the initial 
system Q can be expressed in terms of a certain root of the polynomial S. 

Remark 5. The polynomial S can be found by computing the resultant TZu-^ {h^^^ , r) 
and then dividing it by the polynomial (s^ + l)'^^s"'(a;i3 + yi3sf{xi4 + yi4s)^ 83'^ . 

Remark 6. Since s — tan(p2 — tan 1^93, a complex root of S leads to complex Euler 
angles having no geometric interpretation. Thus, only real roots of S must be 
treated. 

2.8. Structure recovery. Let sq be a real root of S. Then, we can recover the 
matrices i?*^^-', i?'^'^-' and [t]x in closed form as follows. 

We first propose a simple numerically stable algorithm for finding the €,2- and 
Ma-components of the solution. It consists of two steps. First, one finds all real 
roots of the univariate 6th degree polynomials /i^^-'(u2, so) and h^^\u3, so) de- 
fined in (jl8p . Then, the solution (^2,01 ""3,0: sq) corresponds to a minimal value 
of |/i™"^(ti2, U3, so)|j where U2 and U3 run over the obtained roots. 

Now we find the values 

Uj.o = Mj,o/2 - sign(uj,o)^ ("j,o/2)2 4- 1, j = 2, 3, 
subject to \uj^o\ < 1- After that, we obtain Vj^ by and wj^o by 
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where a-"'\ 6p\ are defined in (fT4|) . We also compute the values q, v'^ q 
and Q by After that, we find the entries of R^^^ and [i]x by (O and pU)) 
respectively. 

Let us forget for the moment about the third camera and denote by i?+ — 
R^'^Ku2fi,V2fi,W2fi), ^ R'-'^\u2^q,V2 q,w'2^q), t"^ = ± {t^ ty t^) . Then, as is 
well-known, there are four possible relative positions and orientations for the second 
camera: | t+), | t~), {R~ \ t+) and {R~ \ t~). Let the true configuration 
correspond to (i?^^^ | t). Since the point Pi must be in front of all the cameras, it 
follows that, first, zp^ > and, second, 

{-xo,)Rfj^ + {~yo,)R^23 + {zp, - zo,)R^3^ > 0. 
Denote by ci — t^/R^^, C2 = ciR^^. Then, 

• if ci < and C2 < i^, then i?'^) = R+ ^ t = t+; 

• else if ci > and C2 > t^, then i?^-^' = t = t^; 

• else if tx/Ri3 < and t^R^^/Ri^ < t^, then = R^ , t = t+; 

• else = i?-, t t-. 

Similarly, we find the true relative position and orientation {R^^^ \ at) for the 
third camera. Here a is either +1 or —1. 

After we have found R^^\ R'^^^ and t, the coordinates of O2, O3 and Pi can be 
recovered as follows (recall that the baseline length d = I01O2 fixes the overall scale, 
see Remark [T|): 

d d d 

x'siitx — tzX'21) ZO3 ZO3 

^03 = -TTT , , N ZO2, X03=tx—, yo^=ty—, (35) 

tx — tzX2i ZO2 /oo\ 

zPi = xp^ = xiizp., yp^ = yuzp^. (36) 

Xli X21 f'z 

The true solution to Problem [1] which is assumed to be unique, corresponds to 
a root of S that minimizes the reprojection error 

3 4 

e = ^ ^(a;j» - Xj,)^ + (yj^ - yj^f, (37) 
where the perfectly matched points (xji yji Vj"^ are defined by 



(38) 



As a result, we have obtained a unique solution to Problem[T] which, first, satisfies 
the epipolar constraints ^ and, second, minimizes the reprojection error p7p . We 
must yet multiply the obtained coordinates ([M)) . ([55]) and by pj , where pi 
is defined in ([T|), in order to return to the initial coordinate system. Finally, we 
note that the matrix pjR^^^pj encodes an information on the initial jth camera 
orientation, e.g. the last column of pjR'^^^pj is the z-axis of Ojxyz. 




A SOLUTION TO THE 4P3V POSE PROBLEM IN CASE OF COLLINEAR CAMERAS 11 



' LiHartley ' 



■ New4pt I 



' LiHartley New4pt | 





(a) Generic configuration 



(b) Planar configuration 



Figure 2. Rotational error (in degrees) changing as a function of 
image noise (in pixels) 



3. Experiments on synthetic data 

In this section, we compare our algorithm with the five-point Li-Hartley solver [5] 
on the following two sets of synthetic data: 

(1) generic configuration: all simulated scene points are between the planes 
z = I and z — 2; 

(2) planar configuration: all simulated scene points are on the plane z = 2. 

The third camera center O3 varies randomly between |0i02 and |0i02- The 
baseline length d = 0.3 is the same for both sets of data. The field of view equals 45 
degrees. For each configuration we add the Gaussian image noise with a standard 
deviation of one pixel in a 512 x 512 pixel image. Each experiment is run for 100 
trials. 

The jth rotational error is defined by 



-ij) _ 



arccos 



Tr(i?(j")TpT^O)p^.) 



1 



(39) 



where i?^^-* is the true orientation matrix for the jth camera, Tr(M) is the trace of 
matrix M . The translational error is defined by 



etransl = arCCOS 



(40) 



where t is the true translation vector. 

(2) (3) 

The average values of rotational {— {e).J^ + e).J^)/2) and translational errors are 
reported in Figure [2] and Figure [3] respectively. 



4. Discussion of results 



A non-iterative solution to the four-point three- views pose problem has been pro- 
posed for the case of coUinear cameras. A computation on synthetic data confirms 
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noise noise 

(a) Generic configuration (b) Planar configuration 



Figure 3. Translational error (in degrees) changing as a function 
of image noise (in pixels) 

its correctness and robustness. The new algorithm can be used as a hypothesis- 
generator for RANSAC-like schemes. Its advantages are 

• minimal possible number of scene points needed for reconstruction; 

• uniqueness of the solution; 

• good enough behavior under image noise conditions even in case of planar 
scenes. 

A big number of arithmetic operations needed to derive the polynomial S and 
consequently a big computational error is a weakness of our method in its current 
stage. 
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